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O^' After beginning with a short historical review of the concept of dis- 



placed (coherent) and squeezed states, we discuss previous (often forgot- 



■ ten) work on displaced and squeezed number states. Next, we obtain the 

most general displaced and squeezed number states. We do this in both 
the functional and operator (Fock) formalisms, thereby demonstrating 
the necessary equivalence. We then obtain the time-dependent expecta- 
tion values, uncertainties, wave-functions, and probability densities. In 
conclusion, there is a discussion on the possibility of experimentally ob- 
serving these states. 
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1 Background 



The coherent states were discovered in 1926 by Schrodinger [Q], as an example of how 



his wave functions could mimic classical particles. They consisted of Gaussians in 
an harmonic oscillator potential. The widths of these Gaussians were those of the 
ground-state Gaussian. In modern notation, they were 



It is to be noted that, this being before Born's probability interpretation of 
Schrodinger was concerned by the complex nature of the wave function, wondering if 
it was the real part that was significant. 

Lorentz was the person who had been especially bothered by the lack of classical 
properties of eigenstate wave functions. Schrodinger and Lorentz exchanged many 
letters on this subject, and their correspondence is published 0. 

Soon after this, Kennard published a paper on quantum motion. There he 
described what are, in modern parlance, squeezed states. They follow the classical 
motion, they are Gaussians whose widths are not that of the ground state, and the 
widths and uncertainty products oscillate with time. However, relatively little notice 
was given to this paper in later times 0]. 

In the 1960's the emergence of quantum optics was a fertile background for the 
modern development of coherent states by Glauber, Klauder, and Sudarshan ||5|-0. 
Using boson operator calculus, these states could be defined as either the eigenstates 
of the destruction operator. 



or equivalently as the state obtained by operating on the ground state by the dis- 
placement operator: 




(1) 



a\a) = a\a), 



(2) 




(3) 



D{a) = exp[aa 



a*a] = exp[— |ap/2] exp[aa^] exp[— a*a] , 



(4) 
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the last equality coming from a Baker- Campbell-Hausdorff (BCH) calculation. These 
states are 

CO n 

n=o V n! 

Use of a generating formula for Hermite polynomials shows that this yields the same 
wave function as Eq. (|l]), with the identifications 



V2 



(6) 



Squeezed states were rediscovered and elucidated by a number of people 0-[|T2 
(Ref. |11|] was where the term "squeeze" was invented. See Ref. for a discussion 
of this history.) In operator form, these states are created by 



D{a)S{z)\^) = \a,z) 



(7) 



where the squeeze operator is 
S{z) = exp 



-za'^a^ z aa 
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By BCH relations the squeeze operator can be written as 



z = re^'^ = zi + iz2 ■ 



S{z) = exp 



exp 



^e*^(tanh r)a^a^ 



^e*'^(tanh r)a"fa^ 



/ 1 
V cosh i 



exp 



e *'^(tanh r)aa 



n=0 



(8) 



(9) 



X exp 



--e '"^(tanh r)aa 



The functional forms of D and S are 14 



(10) 



D{a) = exp[—ixopo/2] exp[ipox] exp[— xo5] 



s = s-'/^ 



exp 



iz2 sinh r 
'2^~S 



[X 



exp[— (ln5)(a;9)] exp 



iz2 sinh r 
27 s 



id' 



where 



S = cosh r H — - sinh r = cosh r + cos sinh r = cos^ — he'' sin^ — 
r 2 2 



(11) 
(12) 

(13) 



and one should recall the operator definitions on a function h: 

-- h{x + c) 



exp[cd]h{x) 
exp [r{xd)]h{x) 
exp [c{d'^)]h{x) 



h{xe^) 
1 



POO 


{y-xf 


/ exp 

/ — oo 


4c 



h{y)dy 



(14) 
(15) 
(16) 



[47rc]V2 

Using these functional forms for D and S the most general squeezed wave function 



is 14 



'4^3 



D{a)S{z)% 

1 exp[-z3:opo/2] 
7rV4[5(l+i2fi:)]i/2 



exp 



where 



K = 



2:2 sinh r 
2r5 



Setting z to be real and positive, yields the most commonly studied example: 

{x - Xof 



exp 



tpox 



s = e 



.(17) 



(19) 



2 The Forgotten Displaced Number States 

In the 1950's, before the modern work on coherent states started, there was a short 
flurry of activity by four authors, which has essentially been forgotten. These papers 
were by Senitzky |]T^ (perhaps the best known), Plebanski [|T^, Husimi [T^, and 
Epstein |]T^ . These authors investigated whether there are other wave packets, besides 
the coherent-state packets, which keep their shapes and follow the classical motion. 
They found that there are. In our modern parlance, they found that any displaced 
number state follows the classical motion and keeps its shape. 
In particular, given that the number-state wave functions are 



exp 



X 



HJx) 



7rV4[2"n!]V 



2 ' 



(20) 



the displaced number states are found by applying the functional form of the dis 
placement operator onto the number wave functions: 



exp[-ixoPo/2]exp[ipox] 
ipcs{n) = D{a)ipn = exp 



TT 



[X - Xo) 



[2"n!]i/2 



(21) 



An equivalent result can be obtained in Fock notation by writing 



D(a)\n) = exp[— q;^/2] exp[a;a^] exp[— a*a]|n) 



oo „ fc n 



-p[-«iV2ii:^i:^^ 

fc=0 j=o J- 



{n — j + k)\n\ 



1/2 



(22) 

\n-j + A;). (23) 



_{n - j)\{n - j)\ _ 

Then by changing from the Fock state basis to the wave function basis, interchanging 



the order of the summations, and using the sum 49.4.2 of Ref. |19|, one has 



D{a)\n) 



TT exp[—{x 



a2X — ia\a2\ 



(24) 



2-/2 ^■!(^__^-)! 

Changing the sum over j to a sum over k = n — j, the sum is in the form of the 
seventh equality on page 255 of Ref. . Using this one obtains the same answer as 
Eq. (|1]) above. 

Time-evolution can be shown to be described by |T4 



Xo — > xo(t) = {xq cos t + po sin t) 



The wave packets have the properties 



Po ^ Po(^) = (Po cost - Xo sin t) . (25) 



[Ax(t)f = [Ap(t)]2 = (n+l/2) 



(26) 



{H) = {n+l/2) + \a\\ (27) 

In the later excitement over ordinary coherent and squeezed states, the work of 
the earlier authors [|15]-|]18| was, unfortunately, mainly forgotten. In 1973 one paper 
appeared apUying the operator-formalism displacement operator on the Fock number 



state Then, in the 1980's, papers started appearing studying these states ^2 



24 1 , one of which [22] noted the early work by Senitzky and others [15]-|18 



3 Squeezed Number States 



Plebanski ITEIl, in wave-function form, also looked a little into what we would call 



squeezed number states. By using the functional forms in Eqs. (0) and ([12|) of D 
and S O], we find that the most general squeezed number states are 



D{a)S{z)ilJn 
exp[-ixoPo/2] 

7^1/4 [jrj 1/2 



exp 



[X - Xq) 



-T2 + iv^x 



1 jj f x - Xq 

[2"n!]V2^H ^4 



(2^ 



where 
^2 



5(1 + i2K) = cosh r + e*"^ sinh r, 

/I \ / 1 

-i2K] 



(29) 



52(l + z2/t^ 



52(1 + 4^2) 



z2k 



1 + 



52(l + 4/t2) 



^3 
^4 



1 — z sin sinh r (cosh r + e*"^ sinh r ) 
(cosh r + cos sinh r) (cosh r + e**^ sinh r) 



1 — i2K, ( cosh r + e sin 6 sinh r 



-i2tan-l(2K) 



1 + i2K \ cosh r + e**^ sin sinh r y 

5 ( 1 + ) "'^''^ = [cosh^ r + sinh^ r + 2 cos cosh r sinh r] ^''^ . 



(30) 

(31) 
(32) 



By using properties of the jF's. such as that JF3 is a phase, [T2 + -^2] — 2/^1, and 
T*T\ = one can verify that 



dx 'ip*ss{n)ix) tss{n){x). 



(33) 



Note that when n = 0, the final large square bracket of Eq. ( pH]) is unity and the 
wave function reduces to the normal squeezed states. 

Satyanarayana ^3) defined the problem in Fock notation. He obtained that the 
squeezed number states are 



\a,z,n) = D{a)S{z)\n) 
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eM-W?/'AT.Gmn{z) (^^j L(m)('-™)(|ana'— 10 , (34) 

where the Gmn{z) are given by infinite sums. Later discussions followed. One em- 
phasized number states having the displacement operator and squeeze operator being 
applied separately (rather than in combination) . Another was interested in these 
states from the point of view of the parity operator ||2^. A third still addressed the 
topic from the viewpoint of molecular wave packets p7[ ]. 

Now we will show how to sum the Fock problem in closed form, thereby reproduc- 
ing the result of Eq. ([28|). Start with a number state and apply the squeeze operator. 
If 

(35) 



1/2 



d = -e*'^tanhr, 
2 



then 



S{z)\n) = e'^^" i — — e-'^*'^"|n) 
\ cosh r / 



(36) 



/ 1 X i+ata [2] 
j=0 



\ cosh 'I 



nl 



in-2j)\ 



1/2 



n - 2j) (37) 



]^ xn+l/2 



coshr 



[nl 



,1/2 



1^ (-rf*)^'(coshr)^^' 



k=Q ^■ 



where, in the upper limit of the sums, [n/2] is the greatest integer function. 
With the identification 

' ^ 7rV4[2n-2j+2fc(^_2^- + 2A;)!]i/2' 

and using the known infinite sum 49.4.4 of Ref. |19|, 

4fx^ 



00 k 



-m/2-1/2 



exp 



l + 4v 



Hry 



X 



[1 + 47;] 1/2^ 



(3^ 



(39) 



(40) 



one finds 



where 



S{z)\n) 



exp[-x2j^2/2](n! 



vl/2 



7ri/4jry2 [cosh r (cosh r + e^<^ sinhr)]"/22W2 



(41) 



and 



Hn-2j{xy), 



r = — e '"'^ sinhr(coshr + e**^ sinhr), 



coshr 



1/2 • 



(42) 



(43) 



To do the sum in Eq. (^2]), we start with the n is even case, [?^/2] = n/2. By 
changing the variable of the sum to k = n/2 — j and doing a httle algebra, one finds 



But this sum is of the form of equation 49.4.12 of Ref. 



(44) 



(45) 



{2k)\ ^ (2n)! 

Using this and a fair amount of further algebra one finally obtains for the n being 
even case {n = ue), 



S{z)\n = He) 



7rV4[j-,]iA 



exp 



(^3)"^' 

[2"n!]V2 



(46) 



But this is exactly the form of Eq. (|28|) for the special case of no displacement, 
3^0 = Po = 0. Therefore, now applying the displacement operator trivially yields the 
general result of Eq. (|28D, as it should. 

The n being odd case n = no is obtained in the same manner. The only significant 
difference is that instead of using the sum of Eq. (^Sf ), one uses Eq. 49.4.14 of Ref. 



(2n+ 1)! 



)"(1 + l/t;)"+i/2^2n+i[x(l + \lvY^I\ (47) 



The identical functional form is obtained and so we have demonstrated, as it must 
be, that the functional and operator formalisms yield the identical squeezed number 
states. 



4 Time-Dependent Uncertainties 

Calculating the time-dependent uncertainties is most easily done with the boson- 
operator formalism. This is because we know that 

D''{a)aD{a) = a-a , (48) 
S\z)aS{z) = (coshr)a + e''^(sinh r)a^ (49) 
T\t)aT{t) = ae-^* , (50) 

where T is the time-evolution operator 

T{t) = exp[-i//t] = exp[-i(a^a + l/2)t\ = exp[-i(-9^ + x'^)t/2] . (51) 

With this it is straight-forward to calculate the uncertainties, as a function of time, 
for the squeezed number states. 
To begin. 



{xs{t))n = {n\S^D^T^xTDS\n) = Xq cos t+posmt . (52) 



Similarly, 



{Ps{t))n = pocost - xosint , (53) 
{xl{t))n = (2n + l)[(coshr)^ + (sinhr)^ + 2(coshr)(sinhr)cos(2i-(/))] 

+ Mt))l , (54) 
{pl{t))n = (2n + l)[(coshr)^ + (sinhr)^-2(coshr)(sinhr)cos(2i-0)] 

+ {Ps{t))l ■ (55) 
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Therefore, the uncertainty-product as a function of time is 

[Ax(t)]2[Aj9(t)]2 



{n+l/2f 



1 + 4(cosh r) (sinh r) sin (2t 



l + ^(.^-^)'sin^(2t-0). 



(56) 
(57) 



Eq. (p6D agrees with the results of Ref. [^] for t = and Eq. (^) agrees with the 



standard squeezed-state result when n = [^]. 



5 Time-Evolution of the States 



Taking the functional definition of the time- displacement operator, T, in Eq. (|5TD , 
and using BCH relations it has been shown [jl4| that the harmonic-oscillator time- 
displacement operator can be written as 



T = [cos t] ^^"^ exp 



— tantfa; 
2 ^ 



exp [—(In cost){xd)] exp 



-tant(a^) 



(58) 



Using this result, along with the functional definitions in Eqs. (|T^|I6]), and combining 
it with the squeezed state wave function, ipssin) of Eq. (|28|) , one has 



g— itant x-^/2g— ia;opo/2^'^/2 



TT' 



3/4[isintJ^i2"+in!]i/2 



where 
I -- 



dy exp 



{y — xj cost)' 



\y- x^fT^^iV^y 



y-xp 



i2tant 2 

This integral can be evaluated with the aid of 7.374.8 of Ref. [^, yielding 

1 



2"n! J 



^J-45[A]V2^ 



(59) 



(60) 



exp 



x^ T2 cos t -|- i sin t {xqT2 + ipo) cos t Xq JF2 cos t 

~\~ X 



exp 



2 B 

pI /isint 
" 2" V B 



XoPo 



B 



jFo sin t i 

— h - 

B 2 



5 



(61) 
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where the Ti are given in Eqs. (p9|) - (0), and 



A 
B 

m 



22sint\ /B - i2smt/J^l 

~nB] 



B 



cos t + sin t, 

xo cos t + j9o sin t) . 



X 



(62) 
(63) 
(64) 



Appropriate hmits can be checked. For the case t = 0. one finds ss{n){x-, t = 0) = 
'^ss{n){.x), the function given in Eq. (|28D. Also, for the non-squeezed case of z = or 
jFj = 1, one finds the appropriate time-evolution of displaced number states: 

vi/..(.)(a:,t,^. = l) = ^^T7IJ^;^^n(X(t))exp^ 

exp[z{a; — (xq cost + po sint)/2}(po cost — xq sint)]. (65) 



The probabihty density is given by 



Pss(n){x,t) = ^'*,(„)(x,t)^',,(„)(x,t) 



(66) 



[AA 



*]n/2 



-Hn 



X{t) 



2"(n!)[7r55*]V2jr^ " Vj^45*[A*]i/2^ 
X2(t) 



X(t) 



.^iB[AY/\ 



exp 



(67) 



With the aid of integral 2.20.16.2 of Ref. |g, 



dx e-^^'H4bx)H^{cx) = ^^~y;f {b' + c'-pr/'Pn { , 



(6^ 



where the P„ are the Legendre functions, it can in principle be demonstrated that 
the probability density is properly normalized: 



1 



dx Pss{n){x,t). 



(69) 



However, for small n it is easier to verify the normalization on a case by case basis. 
The first three probability densities are 



Pss{0)ix,t) 



1 



TT 



exp 



BB*J^1 



(70) 
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Pss(o){x,t) is the probability density for the ordinary squeezed state. Pss{i){x,t) and 
Pss(2) {x, t) are n — 1 and n = 2 squeezed number states, having two and three humps, 
respectively, in their wave packets. In Figs. 1-4 we show three-dimensional plots of 
Pss(i) vs. X and t for various n = 1 squeezed states. 

Fig. 1 looks at the case Xq — 8, po — 0, and z — In 2. The picket starts with no 
velocity at maximum displacement. It is at its broadest when t — 0, tt, and 27r for 
X — 8, — 8, and 8. In between the locations where the wave packet is broadest, the 
two humps are closer together, narrower, and higher in amplitude. In Fig 2, we show 
the case again having xq — 8 and po = 0, but this time z — — In 2. This basically 
changes the phase by 7r/2. The narrow, peaked wave packets are when t = 0, tt, and 
27r for x = 8, — 8, and 8. The broader and shorter conditions are in between. 

Fig. 3 is again an example with xq — 8 and po — 0, but now z is imaginary: 
z = (In 2) exp[i7r/2] = iln2. One can see from the t = 2t: case that now we have a 
phase shift such that neither extreme configuration is located at the t = an integral 
number times tt positions. Finally, in Fig. 4, we change to a central position and 
maximum velocity to the right at t = 0: = and po — 8; with z — \n.2. Now the 
broadest wave packets are centered at a; = for t = 0, tt, and 27r. The narrow peaked 
packets are when t — t:/2 and 37r/2. 

6 Discussion 

There is hope that, in the not too distant future, displaced and squeezed number 

states can be observed. This optimism is based on the recent work of Wineland's 
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X^ 



BB*Ti 



X^ 



(71) 



Pss[2){x,t) 



AA* 
2[7r55*]V2jr^ 



exp 




group with trapped Be~^ ions [pTl . 



With laser coohng and a series of laser pulses to entangle the electronic and 
motional states they produced (displaced) even and odd coherent states \TL, 



They also have been able to produce a squeezed (but not displaced) ground state and 
number states If all these techniques can be combined then in principle displaced 



and squeezed number states perhaps can be produced. However, this would be no 
easy feat since there are very complicated heating problems involved. 

When these states are produced, they should mimic the features displayed in our 
Figs. 1 to 4. 
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Figure Captions 

Figure 1. A three-dimensional plot of the n = 1 squeezed-state probability density 
Pss(i) a-s a function of position, x, and time, t, for xq — 8, po — 0, and 2; = In 2. The 
spikes near the maxima of Pss(i) are an artifact of the numerical routine. The true, 
smooth maxima follow the tops of the spikes. 

Figure 2. A three-dimensional plot of the n = 1 squeezed-state probability density 
Pss(i) as a function of position, x, and time, t, for xo — 8, po — 0, and 2; = — In 2. 

Figure 3. A three-dimensional plot of the n = 1 squeezed-state probability 
density Pss{i) as a function of position, x, and time, t, for xq = 8, po = 0, and 
z — (In 2) exp[i7r/2] = iln2. The spikes near the maxima of Pss(i) are an artifact of 
the numerical routine. The true, smooth maxima follow the tops of the spikes. 

Figure 4. A three-dimensional plot of the n = 1 squeezed-state probability density 
Pss(i) as a function of position, x, and time, t, for xq — 0, po — 8, and 2; = ln2. 
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